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TECHNICAL NOTE 2369 


TORSION AND TRANSVERSE BENDING OF CANTILEVER PLATES 
By Eric Reissner and Manuel Stein 


SUMMARY 


The problem of combined bending and torsion of cantilever plates 
of variable thickness, such as might be considered for Bolid thin high- 
speed airplane or missile wings, is considered in this paper. The 
deflections of the plate are assumed to vary linearly across the chord; 
minimization of the potential energy by means of the calculus of vari- 
ations then leads to two ordinary linear differential equations for the 
bending deflections and the twist of the plate. Because the cantilever 
is analyzed as a plate rather than as a beam, the effect of constraint 
against axial warping in torsion is inherently included. The applica- 
tion of this method to specific problems involving static deflection, 
vibration, and buckling of cantilever plates is presented. In the 
static -deflection problems, taper and sweep are considered. 


INTRODUCTION 


For analysis of thin solid wings of small aspect ratio such as 
might be utilized in high-speed airplanes and missiles, beam theory is 
no longer adequate. Wings of this type are more nearly plates than 
beams and should be analyzed by plate theory. However, solutions to 
the partial-differential equation of plate theory axe not readily 
obtained, especially for plates of arbitrary shape and loading. In the 
present paper, therefore, a method is described for obtaining ordinary 
differential equations to replace the partial-differential equation. 

The method employs the minimum-potential-energy principle in conjunction 
with the assumption that the chordwise deflection' shape may be repre- 
sented by terms of a power series. The analysis of the present paper 
is limited to the first two terms of this series. The first term 
represents transverse displacement and the second represents twist. 
Together they permit linear chordwise deflection (the assumption usua ll y 
made in wing design) . Use of the first two terms in the series leads 
to two ordinary differential equations that define the spanwise vari- 
ation of the transverse displacement and rotation. If results of 
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greater accuracy are required, additional, terms in the power series may 
he included (that is, quadratic, cubic, etc.) with a corresponding 
increase in the number of ordinary differential equations obtained. 

Solution of the two ordinary differential equations, subject to 
boundary conditions which arise naturally in the minimization procedure, 
gives the bending deflection and twist at any cross section. The 
stresses may be obtained from the deflections by using the well-known 
equations of plate theory. The order of accuracy of the stresses will 
generally be less than that of the deflections. Because ordinary plate 
theory is employed, the applicability of the analysis is limited to 
plates in which the order of magnitude of the plan- form dimensions is 
greater than ten times the order of magnitude of the thickness. To 
extend the applicability of such an analysis to other dimensions, the 
effects of transverse shear deformation must be included. 

The derivation of the differential equations and boundary condi- 
tions is presented and their application is discussed, and then specific 
problems involving static deflection, vibration, and buckling of canti- 
lever plates are solved. In the static-deflection problems, taper and 
sweep are considered. 

Since completion of the present paper, there has come to the 
attention of the authors a recent paper (reference 1) which presents 
essentially the same ordinary differential equations that are presented 
herein for static -deflection problems. The derivations of these ordinary 
differential equations were done by different methods, and different 
specific problems were solved. 

The basic differential equations presented herein were obtained in 
June 1948 while the fir at -named author was temporarily at the Langley 
Laboratory, and the work was continued by correspondence. 


SYMBOLS 


a, b, i, j parameters specifying taper variation 

c local chord of plate 

h local thickness of plate 

Z length of plate measured perpendicular to root 

m mass per unit area 

A 


P 


lateral load per unit length 
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t 

w 

X, 7 , Z 
D 

E 

P 

T 

w, e 


^" 2 ? a 3 

V n^ 

V P 2 • ’ 
B l> V s 3^ 

c x (x), c 2 (x) 

*1 

Mi 


M, 


M 


M 


x 


M 


x/ 


V 


parameter specify in g twisting -moment distribution 
constant applied twisting moment per unit length, 
transverse deflection, positive in z-direction 
coordinates defined in figure 1 

/ Eh.3 

local flexural stiffness — — 

\l2(l - H 2 ) 

modulus of elasticity of material 
tip shearing force 
tip torque 

functions of x appearing in assumption 
w = V7(x) + y0(x) 

coefficients in differential equations 


functions defining plan form (see fig. l) 

X 

variable obtained by transformation Xi = 1 - by 
tip bending moment 
higher-order moment of stresses 
bending moment per unit width 

exter nally applied tip bending moment per unit width 
twisting moment per unit width 

externally applied tip twisting moment per unit width 
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V 


N x 

normal force per unit width 

N x > N, 
x 0 *1 

constants appearing in normal-force-distribution 

2 y 

equation = IT + — N„ 

x x 0 c T. 

Qx 

shearing force per unit width 


externally applied tip shearing force per unit width 

Iy, Ky 

modified Bessel functions of order V 

j V 

X 

aspect-ratio parameter ( (l - p) J 

P 

Poisson' s ratio 

II 

©- 


CD 

frequency of torsional vibration 

A 

angle of sweep 

U 

energy 

a x 

normal stress 

T T 

xy> xz 

shear stresses 

Subscripts: 


St V 

according to St. Venant torsion theory 

Subscripts 

x and y on w denote partial differentiation with 


respect to x and y, respectively. 

Superscripts : 

h homogeneous solution 

pi particular integral 

Primes denote differentiation with respect to x. 
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AHALTSIS 


The structure considered in the present paper iB a thin, elastic, 
isotropic plate of gradually varying thickness and chord, as shown in 
figure 1. The loading may consist of distributed lateral forces and 
torques, spanwise normal forces acting in the midplane of the plate, 
and tip shears and torques. In addition, the plate may be undergoing 
simple harmonic motion. The potential-energy expression for such a 
plate in its deformed position will now be presented. The aforemen- 
tioned assumption of linear chordwise deformation w ill then be incor- 
porated, and finally the potential energy will be minimized by means of 
the calculus of variations. 

The strain energy of bending is given by the following expression: 


rr _ 1 

^strain 2 J 

pi pc 2 (x) 

D(x,y) ["/ 

2 

'V, + w \ + 

^ xx yyj 



0 c-jjx) 

L 


2(1 

- ( v xy 2 

- v «vj 

J dx dy 

(1) 


where w is the transverse deflection and 


D(x,y) = 


Eh3 

12(l - p 2 ) 


where h is the local thickness which is a function of x and y. 
The potential energy of transverse loads of intensity p is 

pi p c 2 (x) 

Up - -J J p(x,y)w dx dy 
0 c x (x) 

The potential energy of the middle-plane spanwise forces is 


( 2 ) 


% = 


- 1/7 


Cp(x) 


Vx 2 dx dy 


c x (x) 


(3) 



6 


NACA TN 2369 


The potential energy of the tip forces, moments, and torques is 

“tip - I c 5 (-V + "A - Vv) x= j ^ (4) 

If the plate is undergoing simple harmonic motion of circular fre- 
quency cd, and w(x,y) is the deflection shape at the time of maxi mum 
deflection, the potential energy due to inertia loading is 




1 

2 



“A (x,y)cu 2 w 2 dx dy 


( 5 ) 


The total potential energy is defined as the sum of ad 1 the energies 
just listed or 


TI 


total 


= TI 


strain 


+ II 




+ II 


tip 


+ TI 


CD 


If TT fco . fcal were minimized with respect to the deflection w(x,y) 

by means of the calculus of variations, the well-known partial- 
differential equation of plate theory would result as the necessary 
condition the deflection shape must satisfy.. However, if first the 
deflection w is assumed to be of the form 

w = W(x) + y0(x) (6) 

and the potential energy is minimized with respect to W(x) and 0(x), 
two ordinary differential equations are obtained for ¥ and 9. The 
latter procedure is followed herein. It should be noticed that the 
right-hand side of equation (6) is merely the first two terms of a 
power series in y. If greater accuracy is desired, additional terms 

of the series such as y^a(x) and y^7(x) may be used. Substitution of 
expression (6) in the energy formulas (equations (l) to ( 5 )) gives 

“strain " I f [“l"" 2 + 2a 2 W " 0 ” + ^ * 2 < 1 ' * V’ 2 ]** 


( 7 ) 
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H Nx = J (njW' 2 + 2n 2 w , 0' + n 3 0 ,2 ^dx 


( 9 ) 


Utip = -PW( 2) - T0( 2) + M-jW ' ( 2 ) + Mg0*(2) 


( 10 ) 


2 P 2 
_ (a 
TI = -^rr 

CD 2 

V/ 



+ 2s^W0 + s„0 2 ]dx 


e 2 )< 


(11) 


where the coefficients p nJ n^ P, T, and s n are defined 

as follows : 


r c 2 (x) 

v(x) = J t \ 

c x (x) 


D(x,y)y I1_1 dy 


? n (x) * / 


c 0 (x) 


c^x) 


p(x,y)y I1 ~ 1 ay 


n n (x) = J 


c P (x) 

2 N^^dy 
Cl (x) 


f c 2 ( 2 ) _ 

P = J ^dy 

c 1 (2) 


T 


= / , 2 , 1 (v + vK 


c-jlC 2) 


r c p( _ -1 

M n = J . . ^ ^ 


c 1 ( 2 } 


r c 2 (x) 

J n (x) = J , , V x ' y)yI1 

c-jU) 


dy . 
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If "the following variational condition is imposed 

5II total = 6 (^strain + ^ + + ^ip + 3 ^co) = 0 


then 
'Z 

J 0 


m + 


J « ( a l w ")" + ( a 2 0 ")" - Pi _ “ 2 ( s l w + s 2 9 ) + ( n l W ')' + ( n 2 0 ')' 

(a 3 0")" + (a 2 Wy - 2(1 - n)^©')' - p 2 - o> 2 (s 36 + s 2 w) + [nf'y + 


(n #')' 


80 


dx 


(a.-jW")’ + (age")' + n-jW* + n^Bwj - 
Ja.^e n y + ( agW") ' - 2(1 - u)a 1 0' + ^0* + ngW^Js^j + 
jjajW" + a 2 0")8W'J Z + [(a^" + a 2 W")s0^ - 

P 8W( Z) - T 80 (Z) + Mj_ 8W ( Z) + Mg 50 ’ (Z) - 0 ( 12 ) 


At the root (x = 0 ), the following clsmped-edge conditions are 
imposed: 

w(0,y) = v (0,y) = 0 


It follows - that 

W(0) = W’(o) = 0(0) = 0'(O) = 0 (13) 

and, consequently, the variations of these quantities (8W(0), 8W’(0), 

etc.) also vanish. 

At the tip (x = Z), 8W, 80 , GW’, and 80 ' are tahen to he 

arbitrary and, consequently, the tip boundary conditions follow from 
equation (12) in the form 

( a l W ") 1 + ( a 2 0 ") 1 + n l W ' + n 2 0 'l = _P 

— J X—* ~L 


v 
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+ n 3 0' + Hq W 1 = -T (15) 

J x=2 


= -Mq (16) 

x=Z 


(a 3 0" + a 2 W"'j = -Mg (17) 

X - — Z 

The differential equations for W and 0 follow from the vari- 
ational equation (12) in the form 

(aqW")" + (a 20"]" - Pq - a£(sqW + s 2 0) + (nqW')' + (n 2 9 ')' = 0 (±8) 

( a 3 ^")" + ( a 2 W ")" ' 2(1 " ^( a l 0 ')' - P2 " 
co 2 ^ 3 0 + s 2 w) + ^n 3 0‘)' + (ngW 1 ) 1 = 0 (19) 

The problem is to solve equations (l8) and (19) subject to the 
eight conditions given by equations (l 3 ) to (if) • Solution of equa- 
tions (l8) and (19) results in expressions for W(x) and 0(x) . The 
deflection w of the plate is then given by equation (6). 

The stresses may then be calculated by taking appropriate deriv- 
atives of the deflections. The order of accuracy of the stresses will 
be less than that of the deflections since successive derivatives of 
approximate expressions become more and more in error. For this reason 
only those stresses are given which are least subject to error resulting 
from restrictions placed on the deflection function. From known results 
of plate theory 


(a 3 0 n )' + (agW")’ - 2(1 - i^aq©' 


(aqW" + a^") 


tf x (x,y,z) = 


6M X z N x 

h 2 h/2 h 


. . 6M™- z 

t (x,y,z) 

** h 2 h/2 


(20) 


( 21 ) 
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T xz ( x *y* z ) = 


3Qx 

2h 


1 - 


z \2 


h/2 


(22) 


where 


Mx = -D( v xx + ^ w yy) = - D ( w " + Y 0 ") 

Mjjy = (1 - lODw^ = (1 - n)D0‘ 


Q x = ^ - = -D(W" + ye'") - |2( H » + ye") - (i‘- p)^e- 


dx by 


Sx 


dy 


The stresses cr x and are numeric a l l y largest when z = ± — j 

is numerically largest when z = 0. 


SOLUTIONS OF SPECIFIC PROBLEMS 
Outline of Problems Solved 


Solutions are presented for a number of problems involving canti- 
lever plates of various shapes under various loadings. The problems 
may be grouped as follows: 

(A) Rectangular plate of constant thickness 

(1) Tip torque 

(2) Uniform distribution of applied twisting moments 

(3) Torsional vibrations 

( 4) Lateral buckling 
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(B) Symmetrical -plan- form plate with chord and thickness variation 

(1) Symmetric cross section with spanwise variation of chord 
and thickness according to a power law 

(a) Linearly varying chord; tip torque and various span- 
wise distributions of twisting moments 

(b) Chord variations other than linear (solution corre- 
sponding to arbitrary torque loadings left in a formal state 
for a class of chord variations) 

(2) Rectangular cross section with constant chord and expo- 
nential spanwise variation in thickness; tip torque 

(C) Skewed plate of constant thickness and chord under tip loading 
and uniform lateral loading 

The problems of group A were selected because they are simple 
fundamental problems for which solutions obtained by the present method 
can be readily compared with solutions obtained by other methods . Com- 
parison with elementary beam theory is shown for these problems. The 
problems in group B involve tapered plates and are therefore presented 
for their possible application to wing analysis. Problem C was selected 
to show the applicability of the method to swept wings. 

Generally, the plan forms and loadings in the problems chosen are 
those for which the assumption of linear chordwlse deflections might be 
expected to hold - namely, problems involving unswept plates under 
torque loading. A single exception is problem C for which the solution 
presented must be regarded as only a first approximation. 


Rectangular Plate of Constant Thickness 

Foi* a rectangular plate of constant thickness, the flexural stiff- 
ness D is independent of x and y. With the chord of the plate 
denoted by c and with the origin of coordinates at the center of the 
root, the differential equations (l8) and (19) become 


Dc 3 0 rv 
12 


DcW 17 - p 1 - co 2 ( s^W + s 2 e) + (i^W’)' + fae')' = 0 (23) 

2(1 - p)DcS" - p 2 - + s 2 W^ + + /r^W'y = 0 (2^) 



12 


MCA TN 2369 


and the boundary conditions, from equations (13) to (17) > become 

W(0) = W(0) = 0(0 ) = 0 ' (0) = 0 (25) 


f DcW"* + n-jW' + = -P 

> x~ Z 

[^" -2(1 - p)Dc8’ + n 3 0’ + ngW^j ^ = -3 


DcW"( Z) = -M]_ 


Dc- 

12 


■0"( l) = -Mg 


(26) 

(27) 

(28) 
(29) 


For each loading condition the differential equations are solved 
and solutions that satisfy the boundary conditions are obtained in 
closed form. 

Tip torque .- For a plate with a torque T applied at x = 1 , the 
differential equations (23) and (24) become 


DcW 17 = 0 


0 17 - 2(1 - n)Dc0” = 0 


(30) 

(31) 


with the following boundary conditions (from equations (25) to (29)): 

W(0) = W*(0) = 0(0) = 0'(O) = 0 (32) 

DcW" ( z) =0 (33) 


Pc3 

12 


0'"(z) - 2(1 - n)Dc0 1 ( Z) = -T 
DcW" ( Z ) = 0 
— 0 "( 2 ) = 0 


(34) 


(35) 


12 


( 36 ) 
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The differential equations (30) and (31) have the following solutions 
which satisfy the boundary conditions: 


W = 0 


T2 

. , 4Xx 
sinh — 
x 2 

tanh 4X 

. cosh 

2(1 - p)Dc 

2 4X 

kx 1 ' 

1 /. 


(37) 


From this equation it follows that 9 ' is not a constant as in 
the St. Venant torsion theory, in which no constraint against axial 
warping is assumed, but is equal to 


9 ' = — — — — - — (l - cosh + tanh 4X sinh 

2(1 - li)Dcl 2 2 


The twist at the tip (x = 2) is 


0 ( 2 ) = 


T2 


2(1 - |i)Dc ' 


tanh 4X 
4-X 


(38) 


For infinite aspect ratio (X— °°) equation ( 38 ) gives the tip twist 
corresponding to the St. Venant torsion theory 


e(i) 


St V 


T2 

2(1 - (Jt)Dc 


A comparison of the tip twist given by the present theory with the 
tip twist given by the St. Venant torsion theory is presented in 
figure 2 and shows that, for aspect ratios lower than 3 > the tip twists 
given by the present theory are appreciably lower. 

The normal stresses, or so-called bending stresses due to torsion, 
are obtained from equation (20) for the value of 0 given by equa- 
tion (37)* These stresses are 


12Tyz 

1 6 ( 

c^h^ \ 

'l - p\ 


[sinh 


4-Xx 


- tanh 4X cosh 


4-Xx 
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Figure 3 shows the spanwise distribution of the normal stress as esti- 
mated by the present theory for rectangular cantilever plates for values 
of the aspect -ratio parameter X of 1, 2, and 4. This figure shows 
that the normal stress is 0 at the tip (2,y,z) and maximum at the 

extreme fiber at the root This maximum value of the normal 


stress is 



tanh 4X 


Uniform spanwise distribution of applied twisting moments .- For a 
plate with uni form spanwise distribution of applied twisting moments 
(that is, p 2 = t = Constant), the differential equations (23) and (24) 

become 


DcW 1 ^ = 0 


(39) 


e 3 ^ - 2(1 - h)dc 0" = t 
12 


(4o) 


The boundary conditions are the same as those for tip torque, equa- 
tions (32) to (36), except that equation (34) is replaced by 


Dc J 

12 


- 2(1 - ji)Dc0’(Z) = 0 


(41) 


The differential equations (39) and (40) have the following solu- 
tions which satisfy the boundary conditions: 


W = 0 


e = 


t v 


2(1 - n)Dc 


, , 4Xx 
p sinh — t— 

X x^ l 


l 2 V 


4X 


—(tanh 4X + — Ycosh — - 

1+X\ 4X cosh 4 X/\ l 


( 42 ) 
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At x = 2, 


tahh kX _ 1 / 1 

2X 8x2ycosh 4l 

For infinite aspect ratio (A->“) this equation gives the tip twist 
corresponding to the St. Venant torsion theory 



9(1) = 


tl c 


4(l - p.)Dc 


0 ( 2 ) 


t2 2 

St V 4(1 - p)Dc 


A comparison of the tip twist given by the present theory with the 
tip twist given by the St. Venant torsion theory is presented in figure 4 
and shows that, for aspect ratios lower than 5, the tip twists given by 
the present theory are appreciably lower. 

The normal stress at any point and the ma xi mum normal stress are, 
respectively. 


s( W ) - 


sinh 


4\x 


tanh 4X + 


V 


,41x1 

, . , cosh + - — 

4x. cosh 41/ 2 4l 




The spanwise variation of the normal stress as estimated by the present 
theory for rectangular cantilever plates for values of the aspect-ratio 
parameter \ of 1, 2, and 4 is shown in figure 5- 
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Torsional vibrations . - For a plate undergoing torsional vibration 
the differential equation for 0 (equation (24)) becomes 

^e 17 - 2(1 - n)Dc0" - 6=0 (43) 


with equations (32), (36), and (4l) as boundary conditions. The solu- 
tion to the differential equation (43) is 


= A^ sihh ~ + A2 cosh ^ + A3 sin + A^ cos 22 


22 £ 


where 


P 2 . ax2| h + *1 ^ Y* + 1 


l6\ 2 


\“st v/ 


7 2 = 8X 2 


/I + 


it* / m 


l6x 2 ' 


- 1 


ffist vy 


and 


“fet V 


2 6(1 - p.) k 2 D.- 


m.c*Z 


2 7 2 


From the boundary conditions 0(0) = 0 '(O) =0, which are included in 
equation (32), 



a 4 = 


/ 
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From the remaining boundary conditions, equations (36) and ( 4 l), the 
following criterion is obtained: 


1 + 


to qst V 

* CD 


sinh 3 sin 7 + 



32X 2 CDg t v 


CD 


cosh P cos 7=0 


This equation is solved for the fundamental frequency by finding 
the lowest value of the frequency ratio cd^cd^ y that satisfies it for 

a given value of the aspect-ratio parameter X. A comparison of the 
fundamental frequency of torsional vibration given by the present theory 
with that given by the St. Venant torsion theory is presented in 
figure 6 and shows that, for aspect ratios lower than 3, the fundamental 
frequencies given by the present theory are appreciably higher. 

Lateral buck! 1 ng . - For a cantilever plate loaded by a spanwise 
force. 


N Y = N. 


*0 



which is a combination 
the plane of the plate 
are 


of an axial force N Y 

*0 

. The differential 

X 1 


and a bending force in 
equations for this case 


2 ' 

DcW 17 + N cW" + N Y %-0" = 0 ( 44 ) 

x 0 X l 6 


— e 17 - 2(1 - p)Dce" + n x — e " + n x — w" = 0 (45) 

12 x o 12 x i 6 


with the following boundary conditions: 


W(0) = W( 0) = 0(0) = 9' (0) = 0 


(46a) 

c 2 \ 

DcW" + N cW’ + N — 8’) =0 

0 XL 6 Jx=l 

' 

(46b) 
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22^.0 "* - 2(1 - n)Dc0' 
12 


+ N„ 


sle' 


+ N. 


0 12 


*l 6 


c 2 

— W* 


= 0 


x=l 


(46c) 


W "(Z) = 0"( Z) = 0 


(46d) 


Integration of each differential Equation and use of the boundary- 
conditions (46b) and (46c) lead to 


2 

DcW'" + w cW' + 2_0' = 0 (47) 

*0 X 1 6 


3 3 2 

22-0'" _ 2(1 - |J.)Dc0 * + N x — 0* + N x — W* - 0 (48) 

12 12 1 6 


The other boundary- conditions are satisfied by taking 


rr * A, mtx \ 

W = All - cos -gYj 

„ „ nroc\ 

0 = B(1 - cos — yJ 


where n is an integer which represents the number of spanwise buckles. 
Equations (47) and (48) are also satisfied by. these expressions for ¥ 
and 0 if the following stability criterion is satisfied: 


n£ _ Nxo* 2 y n 2 + l6k£ _ NxqZ 2 \ _ l/^V!^ = Q 

, 4 « 2 d A 4 « 2 « 2 d / 3 \ « 2 d J 


(49) 


Equation (49) gives the critical combinations of N and 


N, 


Xq -p 

for a given value of the aBpect-ratlo parameter X. For each value of 
X it is necessary to use the value of n that gives the lowest value 
of K x for a given value of N x > or vice versa. For the present 


0 

problem. 


’V 

n = 1 always gives the lowest values. 
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Timoshenko, in reference 2, has presented a solution for the 
lateral buckling of a strip bent by two equal and opposite eccentri- 
cally applied forces in its plane. As would be expected, the present 
solution differs from that by Timoshenko in that it includes (a) the 
effect of constraint against axial warping and (b) the effect of the 
uniform axial force in reducing the torsional stiffness of the beam. 


Symmetrical Plate with Chord and Thickness Variation 

A class of explicit solutions are presented for torsion of symmet- 
rical cantilever plates with chord and thickness variation. With the 
origin of coordinates at the center of the root the differential equa- 
tions ( 18 ) and ( 19 ) are independent of each other for a symmetrical 
plate in tip torque and distributed twisting moments ( a ^ = 0) . Only 

equation ( 19 ) need be considered; this equation becomes 


(a 3 e-)" - 2(1 - ti)(a 1 0'y - P 2 = 0 


with the boundary conditions 


0(0) = 0 1 (0) = e"(i) = 0 



( 50 ) 


where a^, a^, and P 2 are functions of x defined in the section 

entitled "Analysis." Integration of equation (50) for tip torque alone 
and use of the bracketed boundary condition lead to 

(a 3 e"y - 2(1 - p)a 1 0' = -T (51) 


For applied twisting moments alone, after integration equation (50) 
becomes 



(af")' - 2(1 - p) ai 0' 


(52) 
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Symmetric cross section with algebraic spanwise variation of chord 
and thickness according to a power lav .- Equations (51) and ($2) can b¥ 
solved in closed form when tlie stiffness D (which is proportional to 
■the third power of the thickness) and the chord c vary according to 
the laws 


D - D 0 (l 



c 


- ^ 


tY 


o 

where jSLtSil — } K is a symmetric function of y/c, and Cq 

^(l - p 2 ) 


121 

is the root chord. From the definition of 


*1 


and a„ 


>1 - D o( ] 



a 


3 


-“ol 1 - 



Setting T) 


= ^ leads to 
c 


Do~o^ / rV 2 2 

a 3=— P-T) 12 J. V 2 K( ’>)" d " 
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By use of 



<P = 


ae 

dXi 


P = i + 3d 


q = i + j 


V - k D D 0 


X 0 2 = \ ~^~2 ¥ 1 ~ ^ 

c o 

nl/2 

k D = J , K ^) dT l 

-l/2 



(53) 


m 


equation (51) becomes 


a 

ax^ 




16 Xq 2 „ 12TZ 3 

— P- xf cp = 


D X C V 
^0 0 


(55) 
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and equation (52) "becomes 




12 

D 0 1 c 0 3b 4 





(56) 


It may "be 
cross section, 
table for some 


seen that K = k^ = k^ = 1 for a plate of rectangular 
Values of k^ and k^ are given in the following 
typical cross sections: 


Cross section 

k D 

k X 

Rectangular 

1 

1 

Elliptical 

3* 

32 

2 

Parabolic -arc 

16 

105 

3 

Diamond 

1 

5 


20 



For symmetric sections not given in this table, equations (53) 
and (5*0 may be used. 

This solution is divided into two cases: (a) linearly varying 

chord and (b) chord variations other than linear. 

(a) Linearly varying chordj tip torque and various spanwise 
distributions of twisting moments: In the special case of linearly 

varying chord (j = 1 or p = q + 2) , the solution to the homogeneous 
part of equations (55) and (56) can be expressed in the form 

^ = A l x l Pl + ^l 2 
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where A^ and A^ are arbitrary constants and 


q + 1 / l6\ 0 2 (q + l^ 

Pi + \ / + 


*2 = 


q + 1 / 16 Xq 2 ^ /g + ^ 


The particular solution of equation (55) which applies to tip-torque 
loading is 


cpP 3 - = - 


12T23 


W" 3 1 9 + 


2\ 


"I 


The complete solution to equation (55) therefore is 


cp = cp* 1 + cpP 1 = A-^i^ 1 + AgXp 


32 


12TZ' 


D, 


Ac„3b3| 


0 c 0 b l a - + 2 

b 


V' X1 


Since cp = 

dx 1 


6 = 


Al 


Pl+1 


A 2 


P x + 1 


P 2 + 1 


x/ 2 +1 ♦ A 3 - 


12TZ3 


i6Xo : 


D 0 lc 0 3b3 ( q+ .2 




2 \ 1 - q 
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and A^, 


Ag, and A^ are determined by use of the boundary conditions 


ad) .^-(1) = - b> = 0 

dx-^ dx^ 


The resulting expression for the angle of twist is 



The particular solutions for equation (56) can he found in a 
similar manner when the applied twisting moments pg are known. If 

Pg has the form 



where the value assigned to r defines the distribution of applied 
twisting moments, then 



12p 20 Z 4 

x-jT+l-q (1 - b) r+1 x 1 _q 

+ 1) 

l6\ n 2 i6x n 2 

(r + 2)(r + 1 - q) - q + p 

b 2 0 
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As is done for tip-torque loading, the angle of twist can he found 


to he 


-12p 20 Z j 1 

13 4 ^2 FlI 

D q c 0 h (r + l)\p 2 (l - b) - p^l - b) 


( 7 1 + ^M 1 “ b ) 


7-^r + 1 - q)(l - b) r+1 " q + 7 2 q(l - b) _q 


1 X /l +1 - 1 


P x + 1 


( 7 i + 7 2 ) p i^ 1 " “ 7 i^ r + 1 - q)(i - b ) r+1 q + 


n v P2+1 


r 2 q(l - b) _q 


p 2 + 1 


, r+2-q , 1-q ' 

C 1 ~ 1 X 1 ~ 1 


(57) 


r + 2- q 1-q 


where 


7 1 = 


(r + 2) ( j 


1 6\ n 2 ‘ 


and 


''2 = 


(1 - t ) r+1 



Equation (57) can be used for finding the angle of twist due to 
distributed twisting moment which varies as the rth power of the local 
chord . 
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(b) Chord variations other than linear: The solution to the homo- 

geneous part of equations (55) and (5&) when p ^ q + 2 and" therefore 

j/1 is 


<P 



A l x v 


^tIqCL 

. b 



+ AgK-y 




where I is the modified Bessel function of the first kind of order v 
and Ky is the modified Bessel function of the second kind of order v 
and 


.1 - P._. _ l - i - 33, 

q - p + 2 2(1 - j) 


1 _ q - P + 2 _ t 

a 2 


- J 


The next step is to find the particular integral cpP 1 for the 
torsion load considered from equation (55) or ( 56 ). The complete 

expression for cp is the sum of the homogeneous solution cp* 1 and the 

particular integral cp? 1 . It is then necessary to integrate cp and to 
use the boundary conditions to get the final expression for the twist 9. 

Solution for the twist by means of Bessel functions is straight- 
forward for many values of v. Solutions for cases in which j ^ 1 are 
not carried beyond this point since they involve tabular functions and 
therefore must be worked out separately for any set of values of the 
parameters . 

Rectangular cross section with constant chord and exponential span - 
wise variation in thickness; tip torque .- A case that may be of interest 
is the constant -chord cantilever plate with exponentially decreasing 
stiffness (stiffness is proportional to the third power of the thickness) 


D = D Q e -ax 
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and 


a-^ = DqC6 


-ax 


D 0 c3 

a 3 “ 12 


-ax 


For a plate subjected to tip torque the differential equation (51) is 


( . e") ' 


i6x c 


e-^e* = 


-12T 

D o c3 


( 58 ) 


The solution of equation (58) is 


e 


b-ix box . 3T2 

= A_e + A e d + A + 5 — p 

1 2 3 to 0 c 3 ax 2 


ax 

e 


where 


a + \/a + — 7 t- 


bi = 




2 


For the boundary conditions 0(o) = 0 1 ( 0) = 0"(Z) = 0 , the final 
solution is 


e = 


3T V 


3 2 

to Q c x 


b 0 l b-, Z 

v ~v 


i / boZ aZ\/ bix 

— “t>2 e + ae lie - 1) + 


aZ\/ b2X 


1 ( ^ 1 1 

— (bie - ae lie 


- 1 


1 / ax 

+ t[ e ' 
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Skewed Uniform Plate with Tip and Lateral Loading 

In this section the differential equations for a skewed uniform 
plate under tip and lateral loading are solved, and equations are 
obtained for the twisting and bending deflections of a cantilever plate 
under tip load, tip torque, tip bending moment, and uniform lateral 
load. 


For a skewed plate of chord c (measured parallel to the root as 
shown in fig. 7) and with sweepback angle A 


c 1 (x) = (2 - x) tan A - 


c 2 (x) = ( 2 - x) tan A + 
and for uniform thickness 


a 2 = Dc 


a^ = Dc( 2 - x) tan A 


a^ = Dc 


2 

12 


C + (2 - x) 2 tan 2 A 


Equations (l8) and (19) become 


Dciw^ + j( 2 - x) tan A 0"J 


r = p 


1 . 


( 59 ) 


Dc 


— + (2 - x)^ tan A 
12 


0 " 


+ [(2 - x) tan A ¥”]" - 2(1 - p)0" = 


(60) 
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and "the corresponding "tip boundary conditions from equations (l4) to 
( 17 ) become 



Two integrations of equation (59) with respect to x and use of 
boundary conditions (6l) and ( 63 ) give 


Dc [w" 


+ ( l - x) tan A 




p ( t)dS dn + P(Z 


x) 


\ ( 65 ) 


Integration of equation ( 60 ) with, respect to x and use of 
boundary condition ( 62 ) give 
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If W" from equation ( 65 ) is substituted into equation (66), the fol- 
lowing equation in 9 alone results: 


Dc 



12 


0 ". 


- 2(1 - p)0' 



tan A « 

(2 - x) 

f I dT i + pU - x) - m.^ 


w- 

_ X T] J 

Si 


(67) 


The solution of equation (66) for the case of uniform lateral load 
(pq = pc, P 2 = P C (Z ~ x ) tan A ^ which satisfies the tip boundary condi- 
tion (64) and the root boundary conditions 9(o) = 0'(o) =0 is 


9 = ArYcosh - 1 j + Ag sinh + 


p tan A 


6d(i - p) 


( l - x)3 - i3 


3 l 2 x 


8\ 2 


P tan A 
2Dc(l - p) 




x)‘ 


2 ] 


T 


M-i tan A 

x + x (68) 


2Dc(l - p) 2Dc(l - p) 


where 


A± = -Ag tanh 4X - 


Z 2 

l2Mg 

P tan A 

l6k 2 cosh 4X 

Dc3 

+ 

Dc(l - p) 


P 

p l tan A 

( ± + Jj 

PI' tan A T M 1 tan A 

2D(1 - p) 

V 8k 2 ) 

Dc(l - p) 2Dc(l - p) 2Dc(l - p) 


Equation (68) gives the angle of twist 0 of a skewed uniform cantilever 
plate under tip load P, tip torque T, tip bending moments Mq and Mg, 
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and uniform lateral load p. (if, for example, bending moments are dis- 
tributed symmetrically over the tip of magnitude M, M^_ = M and M 2 = 0.) 

From equation ( 65 ) "with the root boundary conditions W(0) = W'(o) = 0, the 
corresponding bending deflection W is found to be 



- kly? 



+ 3 tan 2 A Z 2 x 2 

+ 2(1 - p ) >. 2 


+ 



T tan A ^ 
2Dc(l - p) 


tan A (l - x ) 9 - tan A 


2k 


[A-^^sinh 



With W and 9 completely determined, the deflection at any point can 
now be found directly from equation (l) and the stresses can be found 
from equations (20) to (22). 


CONCLUDING REMARKS 


A simplified plate theory applicable to thin cantilever plates of 
arbitrary shape and thickness variation and with arbitrary load has been 
presented. The theory, as presented, iB based on the assumption of 
linear deformations in the chordwise direction. Therefore, it is to be 
expected that good results would be obtained for unswept wings in tor- 
sion. For pure bending, the deflections obtained in this manner would 

be off by a factor of as much as 1 - p 2 as a result of the artificial 
restraint against anticlastic curvature. For problems involv ing prin- 
cipally bending, therefore, it may be desired to extend the present 
theory to give more accurate results by considering a more general 
assunption for the deflection w. The expression 


w = W(x) + 0(x)y + a(x)y 2 


which includes the quadratic term <x(x)y 2 in addition to the linear 
terms would give more accurate results. This more general expression 
for w, when used with the energy method, would lead to three linear. 
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fourth-order, differential equations for the quantities W, 
and a complete set of boundary conditions . 


Langley Aeronautical Laboratory 

National Advisory Co mmi ttee for Aeronautics 
Langley Field, Va., March 6 , 1951 
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Figure 2.- Comparison of the tip twist given "by present theory with that 
given hy St. Venant torsion theory for a cantilever plate subjected to 
tip torque. 
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X 


Figure b.- Comparison of tip twist given "by present theory with that 
given hy St. Venant torsion theory for a cantilever plate with a 
uniform distribution of applied twisting moments. 
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l 


Figure 5.- Spanwlse distribution of the normal stress as est ima ted "by the 
present theory for a cantilever plate with uniform distribution of 
applied twisting moments. 
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Figure 6 .- Comparison of the fundamental frequency of torsional vibrations 
given "by the present theory with that given by St. Venant torsion theory 
for a cantilever plate. 




